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Abstract. In this data paper we present the results of an 
extensive 21 cm-line synthesis imaging survey of 43 spiral 
galaxies in the nearby Ursa Major cluster using the West- 
erbork Synthesis Radio Telescope. Detailed kinematic in- 
formation in the form of position-velocity diagrams and 
rotation curves is presented in an atlas together with HI 
channel maps, 21cm continuum maps, global HI profiles, 
radial HI surface density profiles, integrated HI column 
density maps, and HI velocity fields. The relation between 
the corrected global HI linewidth and the rotational veloc- 
ities V max and Vflat as derived from the rotation curves 
is investigated. Inclination angles obtained from the op- 
tical axis ratios are compared to those derived from the 
inclined HI disks and the HI velocity fields. The galax- 
ies were not selected on the basis of their HI content but 
solely on the basis of their cluster membership and incli- 
nation which should be suitable for a kinematic analysis. 
The observed galaxies provide a well-defined, volume lim- 
ited and equidistant sample, useful to investigate in detail 
the statistical properties of the Tully-Fisher relation and 
the dark matter halos around them. 

Key words: Galaxies: fundamental parameters - Galax- 
ies: kinematics and dynamics - Galaxies: spiral - Galaxies: 
structure 



1. Introduction 

The statistical properties of the Tully-Fisher relation 
(TFr), like its scatter and slope, and in more detail the 
characteristics of dark matter halos around galaxies, like 
their core densities and radii, are of great interest to those 
who study galaxy formation scenarios. The relevant ob- 
servables for the TFr are total luminosities and rotational 
velocities while actual rotation curves and luminosity pro- 
files are required to obtain constraints on the density pro- 
files of dark matter halos. However, the available data sets 
from which these observables can be obtained are often 
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suffering from incompleteness, distance uncertainties and 
inhomogeneous observing and data analysis techniques. 

Furthermore, the interpretation of the observables is 
not always unambiguous. For instance, the rotational ve- 
locity of a galaxy is generally inferred from the width of its 
global HI profile and often this measured width is being 
related directly to the dark matter potential (e.g Navarro 
& Steinmetz 200C| ). However, the width of the HI profile 
is a complicated convolution of the 2-dimensional distri- 
bution of HI in a galaxy disk and the shape and extent of 
its rotation curve as sampled by the HI gas. Needless to 
say one has to exercise caution when relating the observed 
scatter and slope in the TFr to the outcome of numerical 
and semi-analytical simulations of galaxy formation. This 
is especially the case when the observed TFr was con- 
structed to serve as an empirical distance estimator and 
the selection criteria and applied corrections to the raw 
observables were optimized to linearize the relation and 
minimize its scatter. 

To overcome some of the observational issues, we initi- 
ated a program to obtain detailed multi-band photometric 
and kinematic information on individual galaxies in a well- 
defined, complete sample. The nearby Ursa Major cluster 
of galaxies provides a particularly suitable sample. In the 
first place, the galaxies in the cluster are all at roughly 
the same distance. Therefore, there is little doubt about 
their relative luminosities, sizes and masses. Furthermore, 
the Ursa Major cluster contains overwhelmingly gas-rich 
systems and the morphological mix of its galaxy members 
is close to that of the lower density field. A detailed dis- 
cussion on the definition of the Ursa Major cluster is given 
by Tully et al. (1996) (Paper I). The characteristics of this 
cluster will be discussed in more detail in Section 2. 

In this paper we present the results of an extensive 
21cm-line synthesis imaging survey of individual galaxies 
in the Ursa Major cluster using the Westerbork Synthesis 
Radio Telescope (WSRT). The HI data are presented in 
the form of an atlas. Multi-band optical and near-infrared 
imaging photometry is presented in Paper I together with 
B-band and K'-band images of all identified cluster mem- 
bers. Forthcoming papers in this series will use these data 



2 



Marc A.W. Verheijen & Renzo Sancisi: The Ursa Major Cluster of Galaxies. IV ; HI synthesis observations 



to investigate the TFr while the HI rotation curves will 
be supplemented in the inner regions with already ob- 
tained high resolution optical rotation curves to derive 
constraints on the structural properties of the dark mat- 
ter halos and the mass-to-light ratios of the stellar popu- 
lations. 

This data paper is organized as follows. Section 2 de- 
scribes the Ursa Major cluster in more detail and con- 
tains morphological and photometric information on the 
galaxies observed with the WSRT. Data acquisition and 
reduction procedures are explained in Section 3. Section 
4 explains how the corrected linewidth can be matched 
to V max and Vfl a t from the rotation curves through an 
appropriate correction for turbulent motion. Inclinations 
are derived from the optical images, HI column density 
maps and HI velocity fields and Section 5 presents a com- 
parison of these inclinations. The layout of the HI atlas 
is described in Section 6 with details about the various 
elements of the atlas pages. The HI properties of the Ursa 
Major cluster galaxies as a sample are presented in Section 
7. Some concluding remarks are given in Section 8. 

2. The sample 

The nearby Ursa Major cluster as defined in Paper I has 
79 identified members. It is located in the Supergalactic 
plane at an angular distance of 38 degrees from the core of 
the Virgo cluster. It has a recession velocity of 950 km s _1 
and a velocity dispersion of only «150 km s _1 . In Paper I 
and in Tully & Verheijen ( |l997|) (Paper II) a distance of 
15.5 Mpc was adopted. However, new HST cepheid dis- 
tances to local TFr calibrators (e.g. Sakai et al. [2000 ) and 
a new correction formalism for internal extinction (Tully 
et al. 1998 ) now place the Ursa Major cluste r gala xies at a 
mean distance of 18.6 Mpc (Tully & Pierce |2000| ). At this 
distance, 1 arcmin corresponds to 5.4 kpc. The morpho- 
logical mix of the cluster members is made up overwhelm- 
ingly by late type systems and only a dozen lenticulars are 
known members. Morphological and photometric proper- 
ties in the optical and near-infrared of individual galaxies 
are described in detail in Paper I. The galaxy distribu- 
tion shows no concentration toward any core and no X-ray 
emitting intra-cluster gas has been detected. 

It should be noted that individual galaxies in the 
nearby Virgo (Warmels [l988a| |l988b| ; Cayatte et al. |1990|) 
A1367 (Dickey & Gavazzi |1991[ ) Her cules (Dickey [1997] ) 
and Coma (Bravo- Alfaro et al. 200C) clusters have also 
been studied in detail and the effect of the dense environ- 
ment on the properties of the HI disks in these clusters 
has been clearly demonstrated. The HI disks in the cores 
of these rich clusters are in general very small and often 
offset from the optical galaxy. On the other hand, the vol- 
ume limited survey of the Hydra cluster (McMahon 1993) 
does not show such an HI deficiency although some inter- 
esting dynamical substructure has been revealed in this 
system. It should be stressed that the Ursa Major cluster 



is markedly different from these more massive and denser 
clusters. In selecting the Ursa Major sample, these envi- 
ronmental effects are carefully avoided as well as fore- and 
background contamination caused by high velocity disper- 
sions and complex dynamical and spatial substructures. 

Since the Ursa Major galaxies are all at the same dis- 
tance, the effects of incompleteness and uncertain relative 
distances are minimized. A complete sample of 62 galaxies 
brighter than M^~— 16.8, i.e. roughly twice the luminos- 
ity of the Small Magellanic Cloud, was constructed and 
nearly all cluster members were observed with the WSRT. 
In this paper, however, only those 49 galaxies which are 
more inclined than 45 degrees, as derived from the opti- 
cal axis ratio, will be considered for a detailed kinematic 
study. 

Table 1 gives a summary of the positional and morpho- 
logical properties of these 49 galaxies while photometries 
are presented in Table 2, based on a 18.6 Mpc distance. 
There are 3 additional galaxies in the tables which do 
not meet the luminosity (*) and inclination ( l ) criteria 
but happened to be in the same WSRT fields as galaxies 
from the complete sample. Of all those 52 galaxies, the HI 
synthesis data of 30 were fully analyzed. Thirteen systems 
were observed and detected but the HI data of these galax- 
ies are presented in an abbreviated form comprising only 
the channel maps, global profiles and position- velocity di- 
agrams. Two of the smaller galaxies were detected in HI 
but they are confused with the HI emission from their 
more massive companions. Finally, there are 7 galaxies in 
the complete sample which have not been observed or de- 
tected because of their low HI content known from single 
dish observations. These are in general SO or Sa systems. 

Table 1 presents the following positional and morpho- 
logical information: 

Column (1) gives the NGC or UGC numbers. 

Columns (2) and (3) provide the equatorial coordinates 

(B1950) derived from the optical images. 

Columns (4) and (5) give the Galactic coordinates. 

Column (6) provides the morphological type. 

Column (7) gives the observed major axis diameter of the 

25 th mag arcsec -2 blue isophote. 

Column (8) contains the position angle of the receding side 
of the galaxy. For galaxies which are not observed or not 
detected in HI, this is the smallest position angle of the 
major axis measured eastward from the north. 
Column (9) contains the observed ellipticity of the optical 
galaxy image. 

Column (10) gives the inclination i opt as derived from the 
observed axis ratio (b/a). See Section 5.1 for further de- 
tails. 

Column (11) gives the adopted inclination angle as derived 
from several methods described in Section 5. 
Column (12) indicates whether a galaxy has a low (LSB) 
or high surface brightness (HSB) according to paper II. 
Columns (13) and (14) provide the galactic extinction in 
the B-band according to Burstein & Heiles (1984) (BH) 
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Table 1. All galaxies in the Ursa Major cluster brighter than M b '*(B)=— 16.8 and more inclined than 45 degrees. 
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Table 2. Photometries of all galaxies in the UMa cluster brighter than M b ' 4 (B)=— 16.8 and more inclined than 45 
degrees. 
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Galaxies with fully analyzed HI data: 
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Galaxies with partially analyzed HI data: 
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Galaxies with confused HI data: 
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and Schlegel et al. (1998) (SFD) as reported by the NASA 
Extragalactic Database. 

Table 2 presents the following photometric informa- 
tion: 

Column (1) gives the NGC or UGC numbers. 
Columns (2)— (5) give the observed total magnitudes in 
the B, R, I and K' passbands from Paper I. 
Column (6) contains the corrected HI line widths at the 
20% level, used to calculate the internal extinction as ex- 
plained below. 

Columns (7)— (10) present the calculated internal extinc- 
tion corrections in thee B, R, I and K' passbands toward 
face-on A^°, calculated according to Tully et al. (1998): 

A*T° = 7a log (a/b) 

where a/b is the observed axis ratio of the galaxy as an in- 
dication of inclination while 7a depends on the luminosity 
and is calculated according to 

7b = 1.57 + 2.75 ( log W^j - 2.5 ) 

7 R = 1.15 + 1.88 ( log W RJ - 2.5 ) 

7l = 0.92 + 1.63 ( log W RI - 2.5 ) 

7 K - = 0.22 + 0.40 ( log W R j - 2.5 ) 

where W^j T is the distance independent HI line width cor- 
rected for instrumental resolution as described in Section 
3.2, corrected for turbulent motion according to Tully & 
Fouque Q1985Q (TFq hereafter) with W t , 20 =22 km/s as 
motivated in Section 4 and corrected for inclination using 
iadopt from Table 1. For dwarf galaxies with W R j < 85 
km/s and for lenticulars with no dust features, the value 
of 7a is set to zero at all passbands. 

Columns (11)— (14) give the total absolute B, R, I and 
K' magnitudes corrected for Galactic and internal extinc- 
tion and a distance modulus of 31.35 corresponding to a 
distance to the Ursa Major cluster of 18.6 Mpc: 



M A ' 1 = m\ ot 



A\ - A\ 



31.35 



where the Galactic extinction A B is taken from SFD as 
listed in Table 1. Extinction corrections in the other pass- 
bands are made according to the Galactic reddening law 
given by Cardelli et al. ( |1989[ ) as summarized by SFD 
under the Landolt filters in their Table 6. It's given by 
A^/A|= 0.62, A\/A^= 0.45 and A^,/A^= 0.08. 
Column (15) gives the diameter of the 25 th mag arcsec -2 
blue isophote corrected for both galactic and internal ex- 
tinction and projection according to TFq: 



Log(D^) = Log(D 25 ) - 0.22 Log(D 25 /d 25 ) + 0.09 A| 

where d 2 5 is the minor axis diameter at the 25 t/l 
mag arcsec -2 blue isophote and Ag is taken from SFD. 



3. Data acquisition and reduction 

The HI data presented in this paper were obtained with 
the Westerbork Synthesis Radio Telescope (WSRT) be- 
tween 1991 and 1996. The integration times varied be- 
tween lx 12 h and 5x 12 h depending on the required signal- 
to-noise. The angular resolution at the center of the cluster 
is 12" x 16" or 1.08x1.44 kpc at the adopted distance of 
18.6 Mpc. The FWHM of the primary beam is 37.4 ar- 
cminutes or 202 kpc. As a result, often more than one 
galaxy was mapped in a single field of view. The ob- 
served bandwidth was either 2.5 or 5 MHz, depending 
on the width of the global profiles. The observations of 
the NGC3992-group and the NGC4111-group required a 
broad frequency band of 5 MHz and at the same time 
also sufficient velocity resolution for the dwarf systems. 
To comply with the correlator restrictions, those two fields 
were observed only in one polarization (XX) which allowed 
for a velocity resolution of 10 km s -1 but resulted in less 
sensitivity. During the earlier measurements an on-line 
Hanning taper was applied but this tapering was aban- 
doned later to obtain the highest possible velocity resolu- 
tion. The various obtained velocity resolutions (dependent 
on the correlator restrictions) were 5, 8, 10, 20 or 33 km/s, 
corresponding to typical rms-noise levels of respectively 
3.1, 1.9, 2.9, 1.6 and 1.0 mJy beam -1 for a single 12 h ob- 
servation at the highest angular resolution. The data of 
NGC4013 were kindly made avaible by R. Botte ma who 
studied this system in great detail (Bottema 1996 and ref- 
erences therein). 

More details on the observational parameters for each 
field are tabulated in the atlas along with the data. What 
follows is a brief description of the reduction procedures. 

The raw UV— data were calibrated, interactively 
flagged and Fast Fourier Transformed (FFT) using 
the NEWSTAR software developed at the NFRA in 
Dwingeloo. The UV points were weighted according to 
the local density of points in the UV plane and a Gaus- 
sian baseline taper was applied with a FWHM of 2293 
(m) which attenuates the longest baseline by 50%. To 
deal with the frequency dependent antenna pattern, five 
antenna patterns were calculated for each data cube at 
a regular frequency separation thoughout the bandpass. 



Pixel sizes of 5 arcsec in RA and 



sm(<5) 



arcsec in declina- 



tion ensure an adequate sampling of the synthesized beam, 
12" x 12'7sin(5). 

After the FFT, the datacubes and antenna patterns 
were further processed using the Groningen Image Pro- 
cessing SYstem (GIPSY). Several channels at the low and 
high velocity end of the bandpass were discarded because 
of their higher noise. As a result, there are 110 or 53 us- 
able channels for a bandpass of 2.5 or 5 MHz respectively, 
except for the N3992 and N4111 fields which had 110 
channels across a 5 MHz bandpass. All datacubes were 
smoothed to lower angular resolutions of 30" x 30" and 
60" x 60". This facilitates the detection of extended low 
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Table 3. A comparison of the widths and integrated fluxes from the present WSRT survey and from the literature. 
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(<0 



10SD 



WSRT 
VLA 

References 



the authors suggest possible confusion with a dwarf companion. 

flagged by the authors as confused with near companion. 

large correction factor (>1.20) applied for primary beam flux attenuation. 

flagged by the authors as possibly interacting. 

no useful single dish profile available due to obvious confusion. 

line width directly measured from the published HI profile. 

the integrated flux as quoted by the author is a factor 2 larger than is quoted by 
any other source. Therefore, half the integrated flux was adopted from this source, 
synthesis observation with the WSRT. 
synthesis observati on w ith the VLA. 
1: Fisher & Tully toll) 8: Schwarz ( |l985| ) 

2: Appleton & Davies ( |1982f) 9: Thuan fc Mar tin (|l98l| ) 
3: Richter & Huchtmeier (p. 99 1]) 
4: Oosterloo & Shostak (|1993 
5: Huchtmeier & Ri chter (n 
6: Schneider e t al. ( |l992j ) 
7: Jore et al. fll99€|) 



10: Magri ([199J) 
11: Van der Burg ( |l987j) 
12: Gottesman et al.j|l984|) 



13: Van Moorsel (1983) 



14: Grewing & Mebold (1975) 



level HI emission and the identification of the 'continuum' 
channels which arc free from line emission. 



3.1. The radio continuum emission 

The channels free from HI emission were averaged and the 
resulting continuum map was subtracted from all channels 
in the data cube. The residuals of the frequency depen- 
dent grating rings were only a minor fraction of the noise 
in the channels containing the line emission. The dirty con- 



tinuum maps were cleaned (Hogbom 1974) down to 0.3cr. 
The clean-components were restored with a gaussian beam 
of similar FWHM as the synthesized beam. When radio 
continuum emission from a galaxy was detected, its con- 
tinuum flux was determined from the cleaned continuum 
map. In cases of no detection, an upper limit for extended 
emission was derived by calculating the rms scatter in the 
flux values obtained by integrating the noise in each of fif- 
teen elliptical areas enclosed by the 25 th mag blue isophote 
and positioned at various emission-free regions in the map. 

3.2. The HI channel maps and the global HI profiles 

At all three spatial resolutions, the regions of HI emission 
were defined by the areas enclosed by the 2(7 contours in 
the 'dirty' 60" resolution maps. Grating rings and noise 
peaks above this level were removed manually. The se- 
lected regions were enlarged by moving their boundaries 
1 armin outwards to account for possible emission in the 
sidelobes. The resulting masks vary from channel to chan- 
nel in shape, size and position due to the rotation of the HI 
disk. These masks defined the regions that were cleaned 
down to 0.3a at all three spatial resolutions. 

The clean-components were restored with a Gaussian 
beam of similar FWHM as the synthesized beam. The 



global HI profiles were derived by determining the primary 
beam corrected flux in each cleaned region. Since the size 
and shape of the clean masks vary as a function of velocity, 
the uncertainty in the flux densities at each velocity in the 
global HI profile varies as well. The noise on the global 
HI profile was determined by projecting each clean mask 
at nine different line-free positions in a channel map and 
integrating over each of them. 

For further analysis, each profile was divided up in 
three equal velocity bins in which the peak fluxes F^ , 

^mid k an d Mngif were determined for the low, middle and 
high velocity bin respectively. These three peak fluxes were 
then used to classify a global profile shape according to: 



Double peaked 

Gaussian 

Distorted 



ipcak 
low 
ipcak 
low 
TTipcak 

low 
ppeak 
01 Mow 
r^pcak 



< F. 



peak 



< F, 



peak 



high 



peak 
mid 
, ppeak 
^ ^mid 
. ppeak 
^ r mid 
peak 



> Fx 



peak 



high 



< F 



peak 



high 



> Fx 



peak 



high 



F 



mid 



F 



peak 



high 



Box y ■ Mow 

In case of a double peaked profile, the peak fluxes on 
both sides were considered separately when calculating 
the 20% and 50% levels. In all other cases, the overal 
peak flux was used. The four velocities V^°, Vf°^, YfA 



and V%lg h corresponding to these 20% and 50% levels 
were determined by linear interpolation between the data 
points, going from the center outward. In the few cases 
of non-monotonically decreasing edges, this procedure 
tends to slightly underestimate the widths. The widths 
are calculated according to 

w _ V 20% \r20% „j w — \r50% \r50% 

W 2 o - V high - V low and W 50 - V high - V low . 
The systemic velocity is calculated according to 
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V sys = 0.25 ( Vj^+V 



20%_lu50%_ 
low 



\t50% 
' v high~ 



-V 



20% 
high 



Because in interferometric measurements some flux 
may be lost due to the missing short baselines, it is useful 
to compare the widths and flux densities from the WSRT 
profiles with those from published single dish observations. 
However, a meaningful comparison requires that the pro- 
file widths are all corrected in the same way for instrumen- 
tal broadening. In general, the widths that are published 
by various authors were corrected for instrumental broad- 
ening using nearly as many different methods. Therefore, 
the published line widths had to be de-corrected first to 
ensure a uniformly applied correction. The de-corrected 
widths and integrated HI fluxes from the literature are 
compiled in columns (7) — (11) of Table 3 along with the 
results from this study in columns (2) — (6). 
Column (1) gives the NGC or UGC numbers. 
Columns (2,3) and (7,8) give the widths of the global pro- 
files at the 20% levels and the formal uncertainties. 
Columns (4) and (9) give the velocity resolutions of the 
observations. 

Columns (5,6) and (10,11) contain the integrated HI fluxes 
derived from the global profiles. 

Column (12) provides the references to the literature 
sources. 

In case the authors suggest that the single dish profile of 
a particular galaxy may be confused and synthesis data 
on that galaxy do exist, these synthesis data are included 
as well and used in the following comparison. However, 
first it will be explained how the observed linewidths are 
corrected for the different instrumental resolutions. 

3.2.1. Correcting W20 for instrumental broadening 

The most widely used method to correct for broadening of 
the global HI profiles due to a finite instru menta l velocity 
resolution was provided by Bottinelli et al. ( 1990| ). For the 
widths at the 20% and 50% levels of the peak flux they 
advocate the following linear relations: 



W 50 „r 



W 20 
W 50 



5W 20 
5W 5Q 



W 20 - 0.55 R 
W 50 - 0.13 R 



where W20 is the observed linewidth and W20.fi is the 
lincwidth corrected for the instrumental velocity resolu- 
tion R in km/s. This empirical prescription is based on 
comparing linewidths obtained at different resolutions. 

However, the correction method applied here deviates 
from Bottinelli et al.'s method and is based on more ana- 
lytic considerations. It is easy to imagine that both edges 
of an intrinsic global profile, when chopped off at their 
peaks and glued together, approximate a Gaussian with 
dispersion ao- The width at the 20% level of this 'true' 
Gaussian is then given by 



A spectral Hanning smoothing was applied to most of the 
WSRT observations presented in this paper. This smooth- 
ing function can also be approximated by a Gaussian with 
a FWHM equal to the instrumental velocity resolution R 
and has a dispersion <tr 



R 



o-r = 



The dispersion er c of the convolved observed Gaussian is 
then given by 



o~ c — y ^0 ~ r w r 

and the 20% line width of this convolved or observed Gaus- 
sian is given by 



W 20 = a c y/8\n(5) 



V8R5) 




R? 
81n(2) 



So, at the 20% level, the intrinsic width W2o,_r is broad- 
ened to W20 by SW given by 



<5W 2n = 




tW81n(5) 



The broadening SW20 does not only depend on the in- 
strumental resolution R but also on the steepness of the 
slopes of the edges of the profile, expressed by erg. Fitting 
Gaussians to the edges of a profile yields o~ c from which 
(Jo can be calculated given the known value of ctr. The 
equation above can be rewritten using a c instead which 
results in 



(R/o-c) 2 
81n(2) 



However, no Gaussians were fitted to the edges of the new 
WSRT profiles. Instead it is assumed that the slopes of 
the edges of the profiles are more or less determined by 
the turbulent motion of the gas with a canonical velocity 
dispersion of <7o=10 km s . This results in 



6W 2 



35.8 ■ 




W 20 ,_r = cr y/ 81n(5) 



and similarly for the 50% level 
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Fig. 1. A comparison of the present WSRT results with 
pre-existing single dish and synthesis data from the liter- 
ature. 




5W 50 = 23.5 



The differences between Bottinelli et al.'s linear 
presciption and our corrections (A<5W = SW Bot - SW our ) 
are only minor and tabulated below for typical instrumen- 
tal resolutions of the WSRT. 



level 


ASW 




5.0 


8.3 


R (km/s) - 
16.5 19.9 


33.1 


20% 
50% 


2.0 
0.2 


2.4 
-0.3 


1.2 0.2 
-3.1 -4.7 


-7.8 

-12.8 



The larger differences occur for the poorest resolutions 
at which only the broadest profiles were observed. Conse- 
quently, the differences are a negligible fraction of the line 
widths. 

Figure 1 shows the comparison of the widths and in- 
tegrated fluxes derived from the new WSRT global pro- 
files and those from the literature. There are no significant 
systematic differences. The unweighted average difference 
in widths is -0.9 ±2.1 km s^ 1 with a rms scatter of 14 
km s _1 . The unweighted average difference in integrated 
flux is 4.7 ± 3.5 percent with a rms scatter of 25 percent. 
It can therefore be concluded that on average the WSRT 
results are in excellent agreement with the results from 
single dish observations. 



3.3. The total HI maps 

As a next step, the total integrated HI maps were con- 
structed from the cleaned datacubes. The clean-masks 
were used to define the regions with HI emission. Out- 
side these regions, the pixels were set to zero and all the 
channels containing a non-zero area were added to build 
up the integrated column density map. This was then cor- 
rected for attenuation by the primary beam. Although the 
advantage of this procedure is a higher signal-to-noise ra- 
tio at a certain pixel in the HI map, the disadvantage is 
that the noise is no longer uniform across the map. As a 
result, the 3er-contour level in an integrated HI map is not 
defined. Signal-to-noise maps have been made, however, 
using the prescription outlined in Appendix A and the 
average pixel value of all pixels with 2.75 < (jt) < 3.25 
was determined. This average value was adopted as the 
'3cr' level for the column density. 



3.4. The radial HI surface density profiles 

The integrated column density maps were used to derive 
the radial HI surface density profiles by azimuthally aver- 
aging in concentric ellipses. The orientations and widths of 
the ellipses were the same as those of the projected tilted 
rings fitted to the HI velocity field (see Section 3.6.1). In 
the case of a warp with overlapping ellipses, the flux in the 
overlapping regions was proportionally assigned to each 
ellipse. The azimuthal averaging was done separately for 
the receding and approaching halves of each tilted ring to 
reveal possible asymmetries. Pixels in the HI map without 
any measured signal were set to zero. Finally, the entire 
radial profile was scaled by the total HI mass as derived 
from the global HI profile. No attempt was made to correct 
the profiles for the effect of beam smearing. 

This method for extracting the surface density profiles 
from integrated HI maps breaks down for nearly edge-on 
systems; the highly inclined annuli with large major axis 
diameters could still pick up some flux along the minor 
axis due to beam smearing. In such cases, Lucy's (1974) 
iterative deprojection scheme as adapted and developed 
by Warmels Qt988b|) might be preferable. 

Due to the complex noise structure of the integrated 
HI map, no attempt was made to estimate the errors on 
the radial HI surface density profiles. 



3. 5. The HI velocity fields 

The HI data cubes were smoothed to velocity resolution of 
«19 km s _1 in order to obtain a good spectral signal-to- 
noisc ratio. HI velocity fields were then constructed by fit- 
ting single Gaussians to the velocity profiles at each pixel. 
Initial estimates for the fits were given by the various mo- 
ments of the profiles determined over the velocity range 
covered by the masks. Only those fits were accepted for 
which 1) the central velocity of the fitted Gaussian lies 
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inside the masked volume, 2) the amplitude is larger than 
five times the rms noise in the profile and 3) the uncer- 
tainty in the central velocity is smaller than i the velocity 
resolution. 

Due to projection effects and beam smearing, the ve- 
locity profiles in highly inclined systems and in the central 
regions of galaxies may deviate strongly from a Gaussian 
shape. The exact shape depends on the spatial and kine- 
matic distribution of the gas within a synthesized beam. 
Fitting single Gaussians to these usually skewed profiles 
results in an underestimate of the rotational velocity at 
that position. As a consequence, the gradients in the ve- 
locity field become shallower. There are several methods 
to correct for the effects of beam smearing. In the present 
cases, however, the signal-to-noise was in general too low 
to allow a useful application of these methods, and, since 
only a small number of systems were recognized as seri- 
ously affected, the HI velocity fields were not corrected for 
the effects of beam smearing. 

3.6. Rotation curves 

The rotation curves were derived in two ways; 1) by fitting 
tilted-rings to the velocity fields (Begeman 1987) and 2) 
by estimating the rotational velocities by eye from the 
position- velocity diagrams. 



3.6.1. using the velocity fields 

The determination of the rotation curves from the velocity 
fields was done in three steps by fitting tilted rings to the 



velocity field (see Begeman 1987, 1989). The widths of the 
rings were set at | of the width of the synthesized beam 
(i.e. 10", 20" or 40"). 

First, the systemic velocity and the dynamical center 
were determined. In this first step the inclination and po- 
sition angles were the same for each ring and kept fixed at 
the values derived from the optical images. The systemic 
velocity, center and rotational velocity were fitted for each 
ring. All the points along the tilted ring were considered 
and weighted uniformly. In general, no significant trend 
as a function of radius could be detected for the systemic 
velocity and center. The adopted values were calculated 
as the average of all rings. 

Second, the systemic velocity and center of rotation 
were kept fixed for each ring while the position angle, in- 
clination and rotational velocity were fitted. All the points 
along the tilted ring were considered but weighted with 
cos(6*) where 8 is the angle in the plane of the galaxy 
measured from the receding side. Hence, points along the 
minor axis have zero weight. While the position angle can 
be determined accurately, the inclination and rotational 
velocity are rather strongly correlated for inclinations be- 
low 60 degrees and above 80 degrees (Begeman p89| ). As 
a result, the fitted inclinations can vary by a large amount 
from one ring to another. However, a possible trend in the 



inclination with radius due to a central bar or a warp can 
be detected. A change in inclination angle often goes to- 
gether with a change in the more accurately determined 
position angle. 

Third, the rotational velocity was fitted again for each 
ring while keeping the systemic velocity, center of rotation, 
inclination and position angle fixed. Again, all the points 
along the tilted ring were considered but weighted with 
cos(#). The fixed values for the inclination and position 
angles were determined in the second step by averaging 
the solutions over all the rings or fixing a clear trend. For 
nearly edge-on galaxies, the inclinations determined in the 
second step were often overruled by higher values based 
on the clear presence of a dust lane (e.g. N4010, N4157, 
N4217) or the very thin distribution of gas in the column 
density maps (e.g U6667). However, uncertainty in the in- 
clinations of nearly edge-on systems does not significantly 
influence the amplitude of the rotational velocity. 

The results of this 3-step procedure were used to con- 
struct a model velocity field. This model was subtracted 
from the actual observed velocity field to yield a map of 
the residual velocities. In some cases (e.g. N3769, N4051, 
N4088) this residual map shows significant systematic 
residuals, indicative of non-circular motion or a bad model 
fit due to a noisy observed velocity field. As a further 
check, the derived rotation curve is projected onto the 
position-velocity maps along the major and minor axis. 

The errors on the inclination and position angles and 
the rotational velocity are formal errors. They do not 
include possible systematic uncertainties due to, for in- 
stance, the beam smearing. 

3.6.2. using the position-velocity diagrams 

It has already been remarked (see Section 3.5) that beam 
smearing affects the determination of the velocity fields, 
especially in the central regions of galaxies and in highly 
inclined disks. As a consequence, the rotation curves de- 
rived from such velocity fields are underestimated as one 
can see from their projection on the XV-maps. In order to 
overcome this problem, the rotation curves were derived 
directly from the major axis XV-maps in a manner simi- 
lar to that used for edge-on systems (cf. Sancisi & Allen 



1979). This was done by two independent human neural 



networks trained to estimate the maximum rotational ve- 
locity from the asymmetric velocity profiles, taking into 
account the instrumental band- and beam-widths and the 
random gas motions. This was done for both the receding 
and approaching side of a galaxy. The rotation curves were 
then deprojected (also accounting for possible warps) by 
using the same position and inclination angles as fixed in 
the third step described in the previous section. In gen- 
eral, the average rotation curves derived from the XV- 
diagrams are in reasonable agreement with those obtained 
by the tilted ring fits. As expected, significant differences 
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can only be noted for galaxies which are highly inclined 
or have a steeply rising rotation curve. 

From the XV-diagrams it is clear that many galaxies 
have kinematic asymmetries in the sense that the rota- 
tion curve often rises more steeply on one side of a galaxy 
than on the other side (e.g. N3877, N3949). The rotation 
curves as derived from the velocity fields and XV-diagrams 
are tabulated in Table 4 for the approaching and reced- 
ing parts separately. The adopted changes in inclination 
and position angles of N3718 and N4138 are motivated in 
the notes on the atlas pages of these galaxies. The uncer- 
tainties quoted in Table 4 are not 1-sigma Gaussian errors 
but rather reflect fiducial velocity ranges, based on the 
position- velocity diagrams. 

4. Matching HI linewidths to V max and Vfl a t 

In this section, it will be investigated how the linewidth 
correction for turbulent motion can be used to match the 
finally corrected global HI linewidths to the actual rota- 
tional velocities measured from the rotation curves. 

After the correction for instrumental resolution, the 
profile widths are generally corrected for broadening due 
to turbulent motions of the HI gas by applying TFq's for- 
mula 



Table 4. Rotation curves derived from velocity fields and 
XV-diagrams. 



wf 



1-2 



2 W x W t , 



( J£i_ 
■(*) 



where the subscript I refers to the widths at the 1=20% or 
the 1=50% level of peak flux. This formula yields a linear 
subtraction of Wt,i if Wi>W c .i and a quadratic subtraction 
if Wi<W c ,i- Values of W t ,i and W c j are different for line 
width corrections at the 20% and 50% levels. The values of 
W c ,i indicate the profile widths where the transition from a 
boxy to a Gaussian shape occurs. The amount by which a 
global profile is broadened due to random motions is given 
by W ti i=2fc;cr where, for a Gaussian velocity dispersion a, 
fc 2 o=l-80 and fe 5Q =1.18. 

The generally adopted values for W c ,i are 
^,20=120 km/s and W Ci 5o=100 km/s. The more impor- 
tant values of Wt,u however, have been subject of some 
debate among various authors. With our new HI synthesis 
data we can give a meaningful contribution to this debate. 



Bottinelli et al. (1983) came up with an empirical ap- 
proach, based on a minimization of the scatter in the TF- 
relation. They assumed an anisotropic velocity dispersion 
of the HI gas of a x =a y =1.5a z and a velocity dispersion 
perpendicular to the plane of <r z = 10 km/s. They deter- 
mined the values of ki by minimizing the scatter in the 
TF-relation and found fc2o=l-89 and £;5o=0.71, indicating 
deviations from a Gaussian distribution (broader wings). 
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he erroneously identifies with Bottinelli et al.'s results. 



16 



Marc A.W. Verheijen & Renzo Sancisi: The Ursa Major Cluster of Galaxies. IV ; HI synthesis observations 



> 
CN 



25 


-25 

25 


-25 



25 

3 

C 

'w -25 

O 

>" 
CN 

25 

* 
-25 



- W ,,20=^8 - 

- • ° * l°° o - 


- W W0 =14 

• • 

* • • 

\ H ° m° 


: i i i i i i i i | h i i : 

W t , 20 =22 , J 
: • * o : 


1 1 1 1 

* .. / • 

« ° »° o ; 



100 200 



100 200 300 



7 W U0^8 ° Q i 

O - 

: .° : 

• * . •; • 

- ' ' ' - 


. i i i i 1 i i b i 1 i i i i _ 

7 W .,ao=14 g ° 1 
: » .° : 

— • " — 

" 1 1 1 1 1 1 1 1 1 1 1 1 


~—\ — 1 — ) — I — — 1 — 1 — \ — \ — | — \ — ] — \ — 

^ W t20 = 32 °c ° i 

r 7 \> • * * : 
* • ■» 

• • 

• _ 

• 


- 1 O 1 _ 

[ W ,.50=15 g ' J 
: • *.° : 

h £~ « * * 

* • 



100 200 100 200 300 

V fIot (kms- 1 ) V flat (kms- 1 ) 

Fig. 2. Comparison of the global profile widths Wr,i, cor- 
rected for instrumental broadening and random motions, 
with 2V ma xSin(i) (Upper panels) and with 2Vfl a tsin(i) 
(lower panels). The left panels consider Wr^o and the 
right panels Wr^q. Different values of the random motion 
parameters Wt,z are used. Open symbol indicate galax- 
ies with declining rotation curves (V max > Vfl at ) and 
filled symbols indicated galaxies without a declining part 
(V max = Vfl a t). See section 4 for further details. 



Finally, Rhee ( |l996a| ) performed the same investiga- 
tion using 28 galaxies, most of them in common with 
Broeils' (1992) sample. Not surprisingly, he found 



Wg%? = 20 ± 2 
= 30 ± 3 



=8±2 
=18 ±3 



similar to Broeils' result. 

Here, with our new and independent dataset, we fol- 
low the same strategy as Broeils and Rhee by investigating 
which values of Wt,i allow an accurate retrieval of V max 
and Vfl a t from the broadened global profile. For this pur- 
pose we will only consider those 22 galaxies in our Ursa 
Major sample that show a flat part in their rotation curves 
(with a significant amount of HI gas) and that are free 
from a major change in inclination angle. Of these 22, 
there are 6 galaxies with V max > Vflat- Note that both 
Broeils and Rhee used Bottinclli et al.'s prescription to 
correct for instrumental broadening which we are forced 
to adopt here to ensure a valid comparison between their 



and our results. We calculated the values of W t n } ax and 
i^flat f or w j 1 j c j 1 th e average values 

AWgf = ^ £ (Wh,i - 2V max sin(i)) 
Aiyff = i £ (Wrj - 2Vfl at sin(i)) 

become zero. This is done for both the entire sample of 
N=22 galaxies and for the subsample of N=16 galaxies 
with V max = Vfl a t. For the entire sample we find 



Wffi = 22 
<!o = 32 



^ = 5 



= 15 



These values are in good agreement with the (rejected) 
results of Broeils and in excellent agreement with the re- 
sults of Rhee. The values of W^f* are larger than the 
values of W t n J ax because the galaxies with V max > Vfl at in 
our sample have considerable amounts of HI gas at their 
peak velocity in the rotation curve. This gas, rotating at 
V max broadens the global profile somewhat further. If we 
consider only the 16 galaxies for which V max = Vfl a t we 
find 



= W t % = 23 



^ t m 5 a o x = ^ t fl 5 1 = 6 



in agreement with the values of W t n } ax we found when 
using all 22 galaxies. 

Our results are illustrated in Figure 2 where we 
show, for each of the 22 galaxies, the deviations 
AW]? ax (upper panels) and AWaf (lower panels) as 
a function of Vfl a t- Galaxies with V max = Vfl a t are 
indicated by filled symbols, galaxies with 
V max > Vfl a t are indicated by open symbols. The 
upper two panels in each block show the results one 
obtains when using Broeil's adopted values of Wt,w = 38 
and W t ,5o = 14 km/s. 

From the upper panels in the upper block it is clear 
that the maximum rotational velocity as derived from the 
corrected global profiles is severely underestimated when 
using the values of Wt,i derived by TFq and adopted by 
Broeils. This systematic underestimation disappears when 
W ti 20 is decreased from 38 to 22 km/s and W t) 5o is de- 
creased from 14 to 5 km/s. The upper two panels in the 
lower block show that if one is interested in the ampli- 
tude of the flat part, which is smaller than the maximum 
rotational velocity for galaxies with a declining rotation 
curve (open symbols), the average offset becomes less sig- 
nificant simply because the open symbols scatter upward. 
In this case, to obtain an average zero offset, we find sim- 
ilar values for W tj / as those adopted by Broeils. However, 
we find the curious situation that the corrected width of 
the global profile systematically overestimates Vfl at for 
galaxies with a declining rotation curve (open symbols) 
and systematically underestimates Vfl a t for galaxies with 
a purely flat rotation curve (filled symbols). 
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From this we can conclude that, in a statistical sense, 
the maximum rotational velocity of a galaxy can be rea- 
sonably well retrieved from the width of the global profile 
when using Wt,20 = 22 or Wt,5o — 5 km/s. The ampli- 
tude of the fiat part can not be retrieved consistently for 
a mixed sample containing galaxies with declining rota- 
tion curves. Note that we have explored only a restricted 
range of rotational velocities: 80-200 km/s. 

Our results also indicate a non-Gaussian dis- 
tribution of random velocities in the sense that 
Wt t w/W t ,50 1-80/1.18. Interpreting W t ,20 and W t , 5 o in 
terms of velocity dispersions it follows that 

c 2 o = W t: 2o/2k 2 o = 6.1 km/s 
cr 5 o = W ti 5o/2k 5 o = 2.1 km/s 

where ^20— 1-80 and ^50=1. 18 for a Gaussian distribution. 
Recall, however, that we advocate a different correction 
for instrumental broadening than Bottinclli ct al.'s scheme 
used by Broeils and Rhee. With our correction method 
for instrumental broadening we find the somewhat smaller 
values of: 



W, 



22 



1,50 



These smaller values of Wt,i allow to retrieve Vfl a t from the 
global profiles of galaxies with purely flat rotation curves 
and V max for galaxies with declining rotation curves. Ap- 
plying our correction method for instrumental resolution 
and the above-mentioned value of Wt 20 — 22 km/s we find 
an rms scatter in AW20 = O.bW^— V max of 6.8 km/s. 



the actual inclination of a galaxy. In the following discus- 
sion we will briefly address those limitations and make an 
intercomparison of i op t, iffli and iyF- 

5-1- iopt from optical axis ratios 

The most widely used formula to infer the inclination an- 
gle from the observed optical axis ratio (b/a)=q was pro- 
vided by Hubble (1926): 



cos 2 (i opt ) 



where qo is the intrinsic thickness of an oblate stellar disk. 
Holmberg (1946) determined an average value of qo=0.20 
which is still commonly used although it is obvious from 
images of edge-on systems that large variations in qo ex- 
ist. For instance, Fouque et al. ( 1990| ) found qo to vary 
from 0.30 to 0.16 for spirals of morphological types Sa to 
Sd respectively and qo=0.42 for galaxies of type Sdm-Im. 
Apart from the debate on the intrinsic thickness, the ob- 
served axis ratio q itself has limited meaning since it is 
often defined at a certain isophote around which q may 
still vary as a function of radius. From images of edge-on 
disks in the Ursa major cluster (see Paper I) it can often 
be observed that the axis ratio keeps increasing outward 
until the faintest isophotes. An extreme example is NGC 
4389, dominated by a narrow bar and surrounded by an 
extended faint halo. The axis ratios presented in Table 1 
were not determined at a fixed isophote but were choosen 
to represent the stellar disk instead of a bulge, lopsided 
structures or a faint halo. 



5. A comparison of inclinations 

Present day instrumentation allows accurate measure- 
ments of the luminosities and global HI profiles of galax- 
ies. In general, the observed scatter in the TF-relation 
is larger than can be explained by the observational un- 
certainties in these measured parameters alone. However, 
the uncertainty in corrections sensitive to inclination con- 
tribute significantly to the observed scatter. For a sample 
of randomly oriented galaxies more inclined than 45 de- 
grees, an uncertainty of 1, 3 or 5 degrees in the inclination 
angle contributes respectively 0.04, 0.12 or 0.19 magni- 
tudes to the scatter due to the uncertainty in line widths 
alone, assuming a slope in the TF-relation of —10. There- 
fore, it is important to determine the inclination angle of 
a galaxy as accurate as possible and this issue deserves 
some special attention. 

From the photometric and HI synthesis data available, 
three independent measurements of the inclination angle 
of a galaxy can in principle be obtained; i op t from the op- 
tical axis ratio, ini from the apparent ellipticity of the HI 
disk, and iyF from fitting tilted rings to the HI velocity 
field. Each of these methods has its own systematic limi- 
tations which are important to recognize when estimating 



5.2. iui from the inclined HI disk 

Apart from the oblate stellar disk, the HI disk can also be 
used to determine the inclination. In general, the HI disk is 
much thinner than the stellar disk and its intrinsic thick- 
ness is of no concern. However, its patchiness, lopsidedness 
and the existence of warps and tidal tails may complicate 
the interpretation of the results from fitting ellipses to a 
certain HI isophote. Here, no correction for the intrinsic 
thickness of the HI layer was applied. However, the rela- 
tively large synthesized beams of imaging arrays at 21cm 
may smear the observed HI disks to a rounder appearance. 
Therefore, a simple correction for beam smearing was ap- 
plied to our measurements and the inclination of the HI 
disk was determined according to 

d 2 - B 2 
cos 2 (i HI ) = ^F-^ 

where Dhi and dm are the observed major and minor 
axis diameters of the inclined HI disk obtained by fitting 
an ellipse to the outer column density levels. 0d and 0d 
are the sizes of the synthesized beam in the direction of 
the major and minor axis of the HI disk. 
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Fig. 3. Intercomparison of the three independently deter- 
mined inclination angles i op t, im and iyF- Only the filled 
symbols (im < 80° when iyF is involved) are considered 
for the unweighted quantitative assessment. 



5.3. ivF from HI velocity fields 

The inclination angle of an HI disk can also be measured 
by fitting tilted-rings to its velocity field (Begeman |1989|) . 
However, the inclination angle and the rotational velocity 
are strongly coupled and reasonable results can only be 
obtained for inclination angles between roughly 50 and 75 
degrees. This procedure requires accurate velocity fields 
with high signal-to-noise ratios as well as many indepen- 
dent points along a ring. The advantage that velocity fields 
offer is the possibility to identify warps and to check the 
kinematic regularity of the HI disk. For instance, the opti- 
cal appearance of a galaxy may look very regular while the 
outer regions of the HI disk may be strongly warped to- 
ward edge-on (e.g N3726). Such a warp would broaden the 
global profile and an inclination correction based on the 
optical axis ratio would lead to an overestimate of the ro- 



tational velocity when dividing the 'warp-broadened' line 
width by sm(i opt ). 

Note that the inclination measurement of a tilted ring 
may be affected by non-circular motions due to spiral 
arms, bars and lopsidedness. 

5.4. the comparison 

For the comparison between the three differently in- 
ferred inclination angles we considered only those 27 
galaxies with fully reduced HI data for which the 
velocity fields and integrated HI maps are available. 
We excluded the interacting galaxies (N3769, N3893, 
U6973) because their outer isophotes (optical and HI) 
are affected by tidal tails. We also excluded galaxies 
with perturbed or inadequately sampled velocity fields 
(N4088, U6969, N4389), galaxies with excessively patchy 
HI maps (N4102) and obviously lopsided galaxies (N4051). 
These eliminations leave us with 19 galaxies that have 
smooth outer isophotes, well filled HI disks and regular 
HI velocity fields. 

Figure 3 presents the comparison between the three 
differently inferred inclination angles using two different 
values for qo. When calculating mean differences and scat- 
ters using iyF, only galaxies with ini<80° are considered 
because kinematic inclinations of highly inclined galaxies 
are systematically underestimated. The error bars on iyF 
are based on the variations in iyF between the various 
fitted rings but are not considered any further here. 

The upper most panel compares iyF with im- No sig- 
nificant offset is found for the 14 galaxies that meet the 
above-mentioned criteria. Assuming that iyF and ini con- 
tribute equally to the scatter of 3.1 degrees implies that 
the inclination angle can be determined with an accuracy 
of 2.2 degrees from either the velocity fields or from the 
inclined HI disk. Note that the correlation turns up for 
iffl>80° due to the systematic underestimation of iyF for 
highly inclined disks. 

Comparing i opt with iyp and ini does show a significant 
offset of roughly 3 degrees when assuming qo=0.20 (middle 
panels). This offset is biggest toward edge-on as would 
be expected in case of an overestimate of the intrinsic 
thickness. Note that there are several galaxies with an 
observed optical axis ratio less than 0.20 which have been 
assigned an inclination angle of 90°. 

This 3° offset disappears when qo=0.09 is used (lower 
panels) and the rms scatter is reduced to only 1.9 degrees 
for i opt versus iyF but is still 4.0 degrees in case of i op t 
versus ini- In the latter case, the scatter is caused by a few 
nearly edge-on systems for which the higher uncertainties 
have no influence on the deprojection of the rotational 
velocities. 

The adopted inclinations and their errors, listed in col- 
umn 1 1 of Table 1 are best estimates based on all the in- 
formation available for a particular galaxy, including the 
morphology of dust lanes if present. For galaxies which 
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Fig. 4. Layout of the HI atlas pages for the 30 galaxies 
with fully reduced data. All the data for these galaxies are 
presented on two facing pages. Results for the 13 galaxies 
with partially reduced data are presented on a single page 
per galaxy and include only the channel maps, the global 
profile and the XV-diagram. The linear scale is 5.4 kpc 
per arcminute. 



lack fully reduced HI synthesis data, the inclination angles 
were inferred from the optical axis ratios using qo=0.09 for 
galaxies of type Sc and later and qo=0.24 for galaxies of 
type Sbc and earlier. The latter value of qo seemed jus- 
tified by the observed axis ratios of the (nearly) edge-on 
systems N4013, N4026 and N4111 of types Sb, SO and SO 
respectively. Unfortunately, there are not enough suitable 
galaxies available to determine qo as a function of mor- 
phology. 



6. The atlas 

The atlas is presented in Appendix B and consists of two 
parts. The first part presents the data for the 30 galaxies 
which have been fully reduced and analyzed. The second 
part gives a less elaborate presentation of the remaining 
13 galaxies which have been only partially reduced. 

The reduction procedures described in Section 3 were 
applied to the data at all three angular resolutions. How- 
ever, to facilitate the intercomparison of various results 
for a galaxy, it was decided to present the data of a par- 
ticular galaxy at the same angular resolution as much as 
possible. The rotation curves are in some cases a combi- 
nation of the rotation curves at various resolutions, the 
inner parts at the highest spatial resolution and the outer 
parts from data of lower resolution but higher signal-to- 
noise ratio. The channel maps are sometimes presented at 
a lower angular resolution than the other data. 

Figure 4 shows the graphical layout for each of the 30 
galaxies in the first part of the atlas. The data for each of 
the 13 galaxies in the second part of the atlas are presented 
on a single page which contains the tables and notes as 
well as the mosaic of channel maps, the global profile and 
the XV-diagram along the major axis. The contents of the 
various tables and panels are described below. 
Tables and Notes - There are three tables presented for 
each galaxy. The upper left table contains information on 
the observations like date, integration time and correlator 
settings. The lower left table presents some of the quan- 
tities derived from the HI data like global profile widths, 
integrated HI flux, systemic velocity etc. The upper right 
table provides the noise and contour levels for the maps 
in the various panels. The contours in the channel maps, 
21cm continuum maps and XV-diagrams are always drawn 
at levels which are related to the rms noise. 

The notes contain information about specific aspects 
of a particular galaxy like optical and HI appearance. 
Channel maps — The channel maps show how the HI 
emission behaves as a function of velocity. From these 
maps it is possible to recognize the presence of warps, 
non-circular motions and Hi-bridges between interacting 
galaxies like N3769/ 11 35+48 and N3893/N3896. The con- 
tours are drawn at levels of -3, -1.5 (dashed), 1.5, 3, 4.5, 
6, 9, 12, 15, ...a. The value of the rms noise level a is given 
in the upper right table. In each panel, the cross indicates 
the adopted dynamical center of the galaxy. The ellipse 
in the upper left panel is centered on this position and 
the position angle is set at the kinematic major axis of 
the inner regions. The ellipticity represents the inclina- 
tion as derived from the optical axis ratio and the major 
axis diameter is equal to D25. The synthesized beam at 
half power is shown in the lower left corners of the left 
panels. The panel at the lower right shows the subtracted, 
'dirty' continuum map. 

Optical image — The optical image of a galaxy was 
scanned from the blue POSS plates. In the upper left cor- 
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ner, the morphological type according to the RC3 is given. 
CCD images of far superior quality can be found in Paper 
I. An example is given for N3726 in which case the CCD 
image is pasted into the POSS image. 

Radio continuum map - The cleaned 21cm radio contin- 
uum map is plotted at the same scale as the optical image. 
The contours are drawn at levels of -4, -2 (dashed), 2, 4, 
8, 16, 32, 64, ...a. The value of the rms noise a is given 
in the upper right table. The ellipse represents the opti- 
cal image as in the upper left panel in the channels maps. 
The synthesized beam is plotted in the lower left corner 
and the small cross indicates the adopted position of the 
center of rotation. 

Global profile - Each point in the global HI profile gives 
the primary-beam corrected, total HI flux density inte- 
grated over a single channel map. The errors are derived 
as explained in Section 3.2. The vertical arrow indicates 
the systemic velocity (lower left table on the first page) as 
derived from the HI velocity field and does not necessarily 
indicate the center of the global HI profile. A significant 
offset from the center could indicate non-circular motions 
or a strong lopsidcdncss. These anomalies can often be 
traced in the individual channel maps. 

Surface density profile - The open and filled symbols in- 
dicate the azimuthally averaged radial surface densities 
for the approaching and receding sides. The solid line fol- 
lows the average value. At the adopted distance of 18.6 
Mpc, 1 arcmin corresponds to 5.4 kpc. The vertical arrow 
indicates R^s- The profile becomes unreliable for highly 
inclined systems since no correction for beam smearing 
was applied. 

XV-diagrams - The position-velocity diagrams are shown 
for two orthogonal cuts through the adopted center of ro- 
tation along the kinematic major (left) and minor (right) 
axes. The position angles of these two axes are printed 
in the upper right corner of each panel. Note that the 
quoted position angles refer to the positive offset axes. 
Consequently the position angle of the major axis refers 
to the receding side which also can be inferred from the 
channel maps. The vertical dashed line indicates the posi- 
tion of the center of rotation. The horizontal dashed line 
indicates the systemic velocity as derived from either the 
velocity field for galaxies with fully reduced data, or from 
the XV-diagram for galaxies in the second part of the at- 
las. The two vertical arrows show where the ellipse with 
major axis diameter D25 intersects the XV-slice. The hor- 
izontal arrows in the left panel show the systemic velocity 
derived from the global profile and ± ^Wrj 
where Wrj is the width of the global profile at the 20% 
level of peak flux, corrected for instrumental broadening 
(see section 3.2) and turbulent motions according to TFq. 

Contours are at levels of -3, -1.5 (dashed), 1.5, 3, 4.5, 
6, 9, 12 , 15, ...g. The value of the rms noise a is given in 
the upper right table. The cross in the lower left corners 
indicates the angular and velocity resolutions. An offset 



of 1 arcmin corresponds to a projected distance of 5.4 kpc 
from the center. 

The crosses give the projected rotational velocities as 
derived from the tilted rings fit. In some cases, at large 
radii, a cross can be found without any underlying sig- 
nal in the XV-diagram. In such cases, the rotational ve- 
locity at that radius is defined by points in the velocity 
field away from the major axis. The open and filled cir- 
cles indicate the projected rotational velocity estimated 
directly from these XV-diagrams. These points must be 
deprojected using the appropriate inclination and posi- 
tion angles to obtain the actual rotation curve for both 
halves of the galaxy. 

Total HI map — All pixels in the total HI map with a pos- 
itive signal have a greyscale value assigned. Because the 
signal-to- noise ratio along a contour is not constant, the 
'3cr-contour' is not defined. Section 3.3 and the appendix 
explain why and how the noise varies across an integrated 
HI map. 

The second contour in the total HI maps corresponds 
to the average value of all pixels with a signal-to-noise 
ratio between 2.75 and 3.25 and this contour can thus be 
considered as a pseudo 3cr-contour. Wherever a contour 
goes through an area with (-j|) < 3, the contour is plotted 
much thinner. Consequently, the lowest contour, plotted 
at the '1.5ct' level, is plotted thin over most of its stretch. 
The various contour levels in atoms cm~ 2 are given in 
the upper right table. The size of the synthesized beam is 
plotted in the lower left corner. The beamwidths are the 
same as in the channel maps unless specified otherwise in 
the note. The small cross indicates the adopted position of 
the center of rotation (1 arcmin corresponds to 5.4 kpc). 
Velocity fields - The greyscales indicate the pixels where a 
radial velocity was measured. Darker greyscales and white 
isovelocity contours indicate the receding side. The thick 
first black contour adjacent to the white ones indicates 
the adopted systemic velocity. In the ideal case of circu- 
lar motion and no noise, this thick contour should be a 
straight line through the center and coinciding with the 
kinematic minor axis of the galaxy. The isovelocity con- 
tours are plotted with constant velocity intervals as given 
by the upper right table. The synthesized beam is plotted 
in the lower left corner. 

The observed velocity field was modeled by fitting 
tilted rings to it. The orientation and rotational velocity 
of each ring were then used to construct the model veloc- 
ity field. The model velocity field is plotted with the same 
orientation and on the same scale as the observed velocity 
field. The isovelocity contours are plotted at the same ve- 
locities in the observed as in the model velocity fields. For 
nearly edge-on systems, the model velocity field is only 
one or two pixels wide in which case no contours could be 
drawn. 

The residual velocity field was made by subtracting the 
model from the observed velocity field. White contours in- 
dicate positive residuals, black contours indicate negative 
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residuals. The contour levels are ...,-15, -10, -5, 5, 10, 15,... 
km s _1 . 

Tilted-ring fits - The three combined panels show the 
results from the tilted-ring fits to the observed velocity 
field. The upper panel shows the inclination angle, the 
middle panel the position angle and the lower panel the 
rotational velocity. 

The crosses with errorbars in the panels for inclination 
and position angle are the results from the second step of 
the fitting procedure as explained in section 3.6.1. The 
dashed lines, mostly coinciding with the solid lines, 



Table 5. Results from the HI synthesis observations. 



Name 


w 20 ± 


W 50 ± 


Res. 


V h cl± 


/Sdv ± 




. ± 


Rjjj 


R lmp 


ylmp 
rot 


± 


shape 


Vmax 


± 


Vflat ± 


£, 










km s 1 








Jy km 


s- 1 


— in 


J y 


I'M 
( ) 


I'M 
( ) 


— km s 


-l 






— km s 


-1 





£0 
i-i 
O 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(ii) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


(18) 


(19) 


(20) 


(21) 


> 


Galaxies with fully reduced HI data: 




































U6399 


188.1 


1.4 


172.5 


2.9 


8.3 


791.5 


0.6 


10.5 


0.3 


<2.5 




1.68 


1.50 


88 


5 


R/F 


88 


5 


88 


5 




U6446 


154.1 


1.0 


131.9 


1.2 


5.0 


644.3 


0.8 


40.6 


0.5 


<7.2 




2.96 


2.93 


80 


11 


F L 


82 


4 


82 


4 




ri.il lb 


286.5 


1.6 


260.6 


1.8 


5.0 


865.6 


0.9 


89.8 


0.8 


49.7 


5.0 


4.24 


6.22 


167 


15 


F/(D) 


162 


9 


162 


9 




N 37b9 


265.3 


6.7 


230.5 


3.6 


8.3 


737.3 


1.8 


62.3 


0.6 


12.1 


2.9 


4.31 


7.10 


113 


11 


F/(D) 


122 


8 


122 


8 




T \ c c n T 


187.5 


1.4 


178.1 


1.9 


5.0 


973.2 


1.2 


11.0 


0.4 


<2.7 




1.64 


1.50 


86 


3 


R L 


86 


3 


86 


3 




N3877 


373.4 


5.0 


344.5 


6.2 


33.2 


895.4 


3.8 


19.5 


0.6 


35.6 


2.4 


2.19 


2.17 


169 


7 


F L 


167 


11 


167 


11 


Pi 


in 3893 


310.9 


1.0 


277.9 


4.1 


5.0 


967.2 


1.0 


69.9 


0.5 


137.4 


2.9 


3.98 


3.88 


148 


19 


F/(D) 


188 


11 


188 


11 


a 


IN 3917 


294.5 


1.9 


279.1 


2.1 


8.3 


964.6 


1.4 


24.9 


0.6 


<7.2 




2.69 


2.83 


137 


8 


F 


135 


3 


135 


3 


N 

o 


JN3949 


286.5 


1.4 


258.3 


1.7 


8.3 


800.2 


1.2 


44.8 


0.4 


134.1 


3.6 


2.62 


1.35 


169 


8 


F L 


164 


7 


164 


7 


m 

£0 

- 




441.9 


2.4 


413.9 


3.2 


33.1 


1052.3 


2.0 


39.3 


0.8 


50.9 


2.5 


3.32 


3.00 


215 


10 


F 


223 


5 


223 


5 


AT 0(170 

IN oy J 2 


281.2 


1 A 

1.4 


260.7 


5.5 


8.3 


852.2 


1.4 


16.6 


0.4 


<5.8 




1.92 


1.67 


134 


5 


R L 


134 


5 






CO* 


U6917 


208.9 


3.2 


189.6 


1.6 


8.3 


910.7 


1.4 


26.2 


0.3 


<4.4 




2.42 


2.00 


111 


7 


R/F 


104 


4 


104 


4 




T T a c\ o o 


166.8 


2.4 


147.1 


4.5 


10.0 


1061.6 


2.2 


10.7 


0.6 


<2.6 




1.29 


1.02 


81 


5 


R L 


81 


5 






H 
V 


U6930 


136.5 


0.5 


122.1 


0.7 


8.3 


777.2 


0.4 


42.7 


0.3 


<5.8 




3.20 


3.17 


108 


7 


R/F 


107 


4 


107 


4 


e Ursa 


N3992 


478.5 


1.4 


461.4 


2.4 


10.0 


1048.2 


1.2 


74.6 


1.5 


30.2 


7.6 


4.75 


6.67 


237 


9 


F/D 


272 


6 


242 


5 


U6940 


59.3 


3.8 


40.6 


7.8 


10.0 


1118.0 


1.7 


2.1 


0.3 


<1.3 




0.61 


0.33 


37 


4 


R 


37 


4 






T TftOft9 


oon q 


ft ft 
0.0 


1 CO A 


Q 7 
O. / 


8.3 


807.4 


3.2 


10.0 


0.3 


13.4 


1.7 


1.38 


1.00 


171 


7 


R L 


171 


7 








N4010 


277.7 


1.0 


264.1 


1.2 


8.3 


901.9 


0.8 


38.2 


0.3 


16.9 


1.6 


3.36 


2.00 


122 


2 


(R)/F L 


128 


9 


128 


9 




U6969 


132.1 


6.4 


123.5 


2.9 


10.0 


1118.5 


2.4 


6.1 


0.5 


<3.8 




0.95 


0.85 


79 


5 


R 


79 


5 






o 


U6973 


367.8 


1.8 


350.4 


1.2 


8.3 


700.5 


1.0 


22.9 


0.2 


127.5 


2.1 


2.21 


1.50 


180 


8 


F/(D) 


173 


10 


173 


10 


Q 


U6983 


188.4 


1.3 


173.0 


1.1 


5.0 


1081.9 


0.8 


38.5 


0.6 


<5.4 




3.07 


3.00 


109 


12 


F 


107 


7 


107 


7 


cT 

CB 


N4051 


255.4 


1.8 


224.6 


1.5 


5.0 


700.3 


1.2 


35.6 


0.8 


26.5 


2.6 


2.89 


2.33 


153 


10 


R/F L 


159 


13 


159 


13 


c+ 
CD 
-! 


N4085 


277.4 


6.6 


255.4 


7.8 


19.8 


745.7 


5.0 


14.6 


0.9 


44.1 


1.3 


1.94 


1.18 


136 


7 


R/F L 


134 


6 


134 


6 


O 


N4088 


371.4 


1.7 


342.1 


1.9 


1Q O 

ly.o 


756.7 


1.2 


102.9 


1.1 


999 3 

zzz.o 


1 Q 

i .y 


A 9^ 


a i n 


1 7 A 


Q 
O 




1 7*} 


1 A 


173 


14 


f Gal 


N4100 


401.8 


2.0 


380.5 


1.8 


19.9 


1074.4 


1.3 


41.6 


0.7 


54.3 


1.7 


3.45 


4.35 


159 


9 


F/D 


195 


7 


164 


13 


N4102 


349.8 


2.0 


322.4 


8.5 


8.3 


846.3 


2.0 


8.0 


0.2 


276.0 


1.5 


1.16 


0.83 


178 


12 


F 


178 


11 


178 


11 


g. 


N4157 


427.6 


2.2 


400.7 


3.1 


19.9 


774.4 


1.8 


107.4 


1.6 


179.6 


2.3 


4.60 


5.67 


185 


14 


F/D 


201 


7 


185 


10 


CD 


N4183 


249.6 


1.2 


232.5 


1.5 


8.3 


930.1 


1.0 


48.9 


0.7 


<5.8 




3.07 


4.02 


113 


11 


F/D L 


115 


6 


109 


4 




N4217 


428.1 


5.1 


395.6 


3.8 


33.2 


1027.0 


3.0 


33.8 


0.7 


115.6 


2.2 


3.19 


3.17 


178 


12 


F/D 


191 


6 


178 


5 


< 


N4389 


184.0 


1.5 


164.9 


1.6 


8.3 


718.4 


1.2 


7.6 


0.2 


23.3 


1.2 


1.30 


1.02 


110 


8 


R 


110 


8 








Galaxies with partially reduced HI data. 


































a 


N3718 


492.8 


1.0 


465.7 


1.0 


33.2 


993.0 


0.8 


140.9 


0.9 


11.4 


0.4 




6.67 


223 


12 


F 


232 


11 


232 


11 


<< 


N3729 


270.8 


1.5 


253.2 


3.9 


33.2 


1059.8 


1.4 


5.5 


0.3 


18.0 


0.9 




1.00 


151 


11 


F 


151 


11 


151 


11 


— 


U6773 


110.4 


2.3 


91.1 


2.2 


8.3 


923.6 


1.6 


5.6 


0.4 


<2.6 






0.67 


45 


5 


R L 


45 


5 






CD 
rjl . 


U6818 


166.9 


2.3 


141.9 


5.7 


8.3 


808.1 


2.1 


13.9 


0.2 


2.4 


1.0 




1.33 


74 


7 


R/(F) L 


73 


5 


73 


5 




U6894 


141.8 


1.1 


132.2 


1.5 


8.3 


848.6 


1.8 


5.8 


0.2 


<2.7 






0.67 


63 


5 


R 


63 


5 






O 

cr 


N3985 


160.2 


3.7 


88.0 


2.4 


8.3 


948.2 


2.0 


15.7 


0.6 


9.7 


1.4 




0.50 


93 


7 


R 


93 


7 






■r. 
CD 
-! 


N4013 


425.0 


0.9 


395.0 


0.8 


33.0 


831.3 


0.6 


41.5 


0.2 


36.3 


0.8 




6.12 


170 


10 


F/D 


195 


3 


177 


6 




U7089 


156.7 


1.7 


97.7 


3.0 


10.0 


770.0 


1.5 


17.0 


0.6 


<3.4 






1.75 


79 


7 


R L 


79 


7 






o 


U7094 


83.7 


1.7 


71.9 


5.5 


10.0 


779.6 


1.6 


2.9 


0.2 


<2.6 






1.00 


35 


6 


R L 


35 


6 









N4117 


289.4 


7.5 


260.3 


5.2 


10.0 


934.0 


1.5 


6.9 


1.1 


3.7 


1.2 










? 










N4138 


331.6 


4.5 


266.0 


7.8 


19.9 


893.8 


3.9 


19.2 


0.7 


16.7 


4.6 




3.55 


150 


21 


F/D 


195 


7 


147 


12 




N4218 


138.0 


5.0 


79.9 


1.9 


8.3 


729.9 


1.7 


7.8 


0.2 


6.3 


0.8 




0.33 


73 


7 


R 


73 


7 








N4220 


438.1 


1.3 


423.3 


3.3 


33.1 


914.2 


1.2 


4.4 


0.3 


<4.9 












? 













Marc A.W. Verheijen & Renzo Sancisi: The Ursa Major Cluster of Galaxies. IV ; HI synthesis observations 



23 



1 r 



^ 



O 



5 0\r 



cn 
O 



O 

- 1 -2 tr 



1 I 1 I 1 I 1 I 1 4- 1 I 1 I 1 I 1 I 1 I 1 I M 



9, o 

O OD 

6ocf<A + c 



i i i i i i i i i 



1 r :** • 



I I I I I I I I I 



„ <s> ° <P 

°fO„ 



1 r »J-»o- ° 



8S ; 



1 1 1 1 1 1 1 1 1 1 1 1 1 



1 1 1 1 1 1 1 1 1 1 1 1 1 



* . • : e 



o o 



I . I . I . I ' T i I i I i I i I i I i I i H 



-24 -22 -20 -18 

M K , 



SO Sa Sb Sc Sd Sm 

Galaxy Type 



Fig. 5. Correlations between HI mass-to- light ratios and 
absolute if'-band magnitudes (left column) and morpho- 
logical type (right column). Solid symbols indicate HSB 
galaxies and open symbols denote galaxies of the LSB 
type. 



in these upper two panels indicate the final values of the 
inclination and position angles kept fixed when the rota- 
tional velocity was fitted. The resulting rotation curve is 
shown by crosses with errorbars in the lower panel. The 
errorbars indicate the formal errors, as given by the least 
squares minimization algorithm. 

The horizontal arrows in the upper two panels indi- 
cate the inclination and position angles as derived from 
the optical isophotes in the outer regions. The diamonds 
indicate the inclination and position angles as determined 
from the total HI maps. When the total HI maps are very 
patchy, these diamonds are very uncertain. The horizontal 
arrow in the lower panel indicates the rotational velocity 
as derived from the width of the global HI profile corrected 
for instrumental broadening, turbulent motion and incli- 
nation. The adopted inclination is representative for the 
outer parts. The vertical arrow in the lower panel indicates 

Rb,i 
25' 



The solid lines in the upper and middle panels show 
the inclination and position angles that were adopted to 
deproject the radial velocities determined from the XV- 
diagrams. This deprojection results in the rotation curves 
plotted as open and filled circles in the lower panel (same 
symbols as in the XV-diagrams) . Note that although the 
rotational velocities at a certain radius may be different 
for the approaching and receding sides, both sides were 
assumed to have the same inclination and position angles 
at that radius. The solid line in the lower panel shows 
the mean rotation curve derived from the XV-diagram. 1 
arcmin on the horizontal axis corresponds to 5.4 kpc. 
Tabulated data — The various parameters derived from 
the HI data presented here are summarized in Table 5. 
Column (1) gives the NGC or UGC numbers. 
Columns (2)-(5) give the uncorrected widths with formal 
errors of the global profiles at 20% and 50% of the peak 
flux. 

Column (6) gives the instrumental velocity resolution at 
which the global profiles were observed. 
Columns (7) and 8 contain the heliocentric systemic 
velocities and their uncertainties as derived from the 
global profiles. 

Columns (9) and (10) provide the integrated HI flux and 
the uncertainty in Jy km s . 

Columns (11) and (12) contain the 21cm continuum flux 
density and its uncertainty in mJy. In case no continuum 
flux was detected, a 3a upper limit for extended emission 
is given. 

Column (13) gives the radius of the HI disk, Rhi in 
arcmin, at the azimuthally averaged surface density of 
1 M Q pc -2 , measured from the radial surface density 
profiles. 

Column (14) gives the radius R lm P of the last measured 
point of the rotation curve in arcmin. The differences 
between R lmp and Rhi depend on the sensitivity of the 
measurement and the distribution of the HI gas along the 
kinematic major axis. 

Columns (15) and (16) give the rotational velocity of the 
last measured point V lmp and its uncertainty. 
Column (17) contains information on the overal shape 
of the rotation curve; R: rising rotation curve, F: the 
rotation curve shows a flat part, D: the rotation curve 
shows a declining part, L: lopsided. 

Columns (18) and (19) give the maximum observed 
rotational velocity V max and its uncertainty. For galaxies 
with a rising rotation curve (R) V max = V lmp . 
Columns (20) and (21) give the average rotational 
velocity of the flat part of the rotation curve Vfl a t and its 
uncertainty. For galaxies with a flat rotation curve (F) 
Vfl a t = V max may deviate from V lmp because Vfl at was 
averaged over the flat part of the rotation curve while 
V lmp was measured at a single point. 
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Fig. 6. Ratio of Hi-to-optical diameter as a function of K' magnitudes, morphological type and disk scale length. 
The HI diameters were measured at the 1 M Q pc -2 isodensity contour. Solid symbols indicate HSB galaxies and open 
symbols denote galaxies of the LSB type. 



7. HI properties of spirals 



7.1. Global parameters 



The HI survey of the Ursa Major cluster presented here 
provides not only the kinematical information necessary 
for the study of the Tully-Fisher relation and of the dark 
and luminous matter for a well defined sample of galaxies. 
It also serves to investigate the general HI properties of 
disks and to make a comparison with galaxies in the field 
and with galaxies in denser environments. The HI studies 
of the Virgo cluster galaxies by Warmels (1988a, 



and especially by Cayatte et al. (199C) have shown that 
the spiral galaxies in the central parts of the cluster have 
smaller HI disks of lower surface density. In the Hydra 
Cluster McMahon ( 1993Q did not find any such significant 
HI deficiency. She did find, however, a surprisingly large 
number of isolated Hl-rich dwarf galaxies near the center 
of the cluster. Dickey (1997) surveyed the more distant 
Hercules SuperCluster and found a similar HI deficiency 
of spirals near the X-ray gas as in the case of Virgo. 

The Ursa Major cluster differs from those just men- 
tioned. It has no central concentration, no X-ray emit- 
ting gas and contains mainly spirals of late morphological 
types. In many respects its conditions are very similar to 
those of a field environment. For this reason it is useful 
to compare the properties of the Ursa Major spirals not 
only with those of galaxies in dense cluster environments 
but also with those of field gal axies as found in various 
recent studies (se e e.g. Br oeils 1992 , Puche & Carignan 



1991 , Rhee |l996a| , Swaters |l99g) 



Here we give only a brief description of the global pa- 
rameters and of the main properties of the HI disks of 
the Ursa Major galaxies. A more detailed discussion and 
a comparison with results from previous work is beyond 
the scope if this data paper. 



Integral properties and global parameters of spiral galax- 
ies have been derived for a large number of objects from 
single-dish observations (cf. Roberts and Haynes 1994) 



In recent years also synthesis observations (Broeils 1992 



Rhee 1996a) have been used to obtain similar information 
for smaller samples of galaxies. 

The Mhi/L ratios obtained for the galaxies of the Ursa 
Major sample listed in Table 5 are shown here in Figure 5 
as a function of absolute magnitude and of morphological 
type. It is well known (see refs. above) that the Mhi/L 
ratio of galaxies depends on luminosity and morphological 
type. The present sample of galaxies shows a clear increase 
of the HI mass fraction with decreasing luminosity and 
from early to late morphological types. The correlation is 
clearly stronger for the i^'-band magnitudes which is a 
better tracer of the stellar mass. 

7.2. Sizes of HI disks and radial surface density profiles 

Detailed information on the sizes and radial distributions 
of HI disks has been obtained recently from synthesis ob- 
servations of limited samples of field and cluster galaxies 
(Broei ls a nd Van Woerden 1994 , Cayatte et al. 1994 , Rhee 



1996a and 19961 ) . Here we present only some of the main 



results on the comparison of HI and optical diameters, on 
the relation between HI mass and diameter and on the 
radial density profiles for the Ursa Major sample. 

Figure 6 shows the ratio of the HI diameter Dhi (de- 
fined at an HI surface density of 1 M0pc~ 2 ) to the optical 
diameter D 2 5 as a function of luminosity, morphological 
type and disk scale-length. The diagrams do not indicate 
any clear trend or dependence of the diameter ratio on 
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Fig. 7. Correlations between HI mass and the isophotal 
diameters of the stellar (triangles) and HI disks (circles). 
Solid lines indicate fits to the plotted data points w hile th e 
dashed lines represent the fits found by Broeils ( 1992Q . 
Filled symbols indicate HSB galaxies and open symbols 
denote galaxies of the LSB type. Triangles are offset by 
0.3 dex to the left. 



any of those quantities. The spread is large. There may be 
a hint of a slight increase of the ratio from early to later 
types and from more luminous to less luminous systems. 
For almost all galaxies Dhi is larger than D 2 'g. 

As shown in previous investigations (see refs. above), 
there is a tight correlation between HI mass and HI di- 
ameter as illustrated in Figure 7. This implies a nearly 
constant mean HI surface density regardless of size. The 
HI mass correlates also with the optical diameter, but, as 
in previous work, with a much larger scatter. 
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7.3. Warps, asymmetries and interactions 

The radial distributions of the HI surface densities are 
shown in Figure 8. Only galaxies with fully reduced data, 
more inclined than 80 degrees and with Rhi < 1 arcmin 
are considered in order to avoid the most severe cases of 
beam smearing. There is clearly a considerable diversity 
of shapes and intensities. The upper row shows the pro- 
files grouped for galaxies of similar morphological types. 
The dotted lines represent low surface brightness galax- 
ies. No obvious trend with morphological type or surface 



brightness can be discerned. However, in the lower row, 
the profiles are grouped according to the galaxy proper- 
ties as listed in Table 6. In this clear trend is visible 
in the sense that galaxies with high HI surface densities 
in their inner regions are either involved in interactions, 
are lopsided or display a warped HI disk. Especially in 
the case of interacting and strongly lopsided systems, no 
longlived stable gas orbits can be expected and evidently, 
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the cold gas becomes concentrated toward the inner re- 
gions of the disk. Note that some galaxies appear in more 
than one panel. 

Warps are thought to be a quite common feature of the 
outer HI layers of spiral galaxies. But, in spite of attempts 
made in recent years (Bosma, 1991), a good statistics on 
their occurrence does not exist yet. The information on HI 
warps provided by the present survey of the Ursa Major 
sample is somewhat limited mainly because of the already 
noted small radial extent of the HI layers with respect 
to the optical. There are pronounced HI warps like the 
well-known ones found in NGC 3718 (Schwarz 1985) and 
NGC 4013 (Bottema 1995) and those seen in other galax- 
ies (N3726, N3985 N4010, N4157, N4183). In most cases 
the warps are visible in the outer parts beyond the optical 
bright disk and are present in normal, regular, not inter- 
acting disks. But there are also systems, like NGC 4088, 
which are strongly distorted in their optical appearance 
and also in their kinematics. Similar distortions are found 
in clearly interacting systems like NGC 3769. In the sam- 
ple of 43 galaxies (Table 6) there are at least 13 objects 
with clear indications of warping. 

Also asymmetries are thought to occur frequently in 
field spirals (Richter and Sancisi, 1994). In the present 
sample we see a large number of objects (at least half, see 
Table 6) with a lopsided HI distribution and/or kinemat- 
ics. The majority shows kinematical asymmetries: on one 
side of the disk the rotation curve rises more slowly and 
reaches its fiat part at a larger radius than on the other 
side (see for example NGC 3877, 3949, 4051). Note that 
this occurs in the inner parts of non-interacting, regular, 
normal systems. 

Finally, there are four galaxies in this sample of 43 
which have close companions and show clear HI signs of 
tidal interactions. A few more objects (examples NGC 
3718, 4088) show distortions or peculiar structures in their 
density and velocity maps. 

8. Concluding remarks 

This data paper has presented the results of an extensive 
HI synthesis imaging survey with the WSRT of a well- 
defined complete equidistant sample of spiral galaxies in 
the nearby Ursa Major cluster. Figures B.l and B.2 show a 
compilation of all the available HI maps. Individual galax- 
ies are at their proper position on the sky but they are 
individually four times enlarged. Some galaxies had to be 
shifted to avoid overlapping maps. Rotation curves have 
been derived for most galaxies as well as detailed infor- 
mation on the kinematical state of the galaxy disks as 
indicated by the presence of global perturbations, warps, 
interactions and lopsidedncss. Since the galaxies were not 
selected on the basis of their HI size or content, the qual- 
ity of the kinematical data varies widely from galaxy to 
galaxy. Nevertheless, these data will be useful for an analy- 
sis of the statistical properties of the Tully-Fisher relation 



and the rotation curves may be decomposed into contribu- 
tions from the main dynamical constituents like the stel- 
lar and gaseous disks, the bulge and the dark matter halo. 
These issues will be addressed in forthcoming papers. 

Those who would like to use the rotation curves for 
their own purposes are advised not to take the rotation 
curves at face value but to check their validity against the 
actual data and to take notice of the comments provided 
on the atlas pages. It should be realized that the rotation 
curves presented here are derived from the kinematics of 
the HI gas which is assumed to be a tracer of a galaxy's 
potential. There are some caveats to be aware of like the 
projection effects of streaming motions and partially filled 
HI disks of edge-on systems. 

Finally, optical long-slit spectroscopy is already avail- 
able for nearly all galaxies in the sample. These high res- 
olution optical rotation curves will be used to supplement 
the HI rotation curves in the inner regions, allowing for 
better decompositions and maximum-disk constraints. 
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Appendix A: The noise in an integrated HI map. 

This appendix explains how the noise in a total HI map 
can be calculated. A total HI map is usually constructed 
from a 3 dimensional datacube containing a number of so 
called channelmaps. Each channelmap shows an HI image 
of the galaxy at a certain velocity. A total HI map is made 
by adding those channelmaps which contain the HI signal. 
Before adding the channel maps the signal in each chan- 
nelmap should be isolated. When the signal is not isolated 
one merely adds noise to the total HI map because the lo- 
cation of the signal in a channel map varies with velocity 
due to the galactic rotation. The signals can be isolated 
interactively by blotting away the surrounding noise or in 
a more objective way by taking a certain contour level in 
the smoothed maps as a mask. As a consequence of adding 
channel maps with isolated regions, the noise in the total 
HI map is not constant but varies from pixel to pixel. The 
noise at a certain pixel in the total HI map depends on 
the number N of non-blank pixels at the same position in 
the individual channel maps that were added. 

In case the data cube was obtained with an uniform 
taper during the observation, the noise a u in two chan- 
nelmaps will be independent. The noise equivalent band- 
width B" in a uniform tapered spectrum is equal to the 
channel separation b. When adding N uniform tapered 
channelmaps at a certain pixel the noise at the same 
pixel position in the total HI map will be increased by a 
factor >/~N : 

<j u N = VN a u 

Usually, the observations are made using a hanning ta- 
per in which case the noise in two adjacent channelmaps 
is no longer independent. A hanning taper effectively 
smooths the data in velocity by convolving the velocity 
profile at each pixel. If Uj is the pixel value in the \ th uni- 
form tapered channel map, the value Hj in the i th hanning 
tapered channelmap is given by 

Hi = jUi-i + \Ui + \U{+i 

Since the cr"'s are independent and all equal to <r u , 
the noise a 1 } in the i th hanning tapered channelmap can 
be calculated according to 

„h |Yl_u ^2 | / 1 „u\2 | (l„u \2~\ h 




= ^ g u = 0.61 a 



In this case the noise equivalent bandwidth B h for a han- 
ning tapered spectrum is given by 

B h = i£ B u = 2 6? B u 



As a consequence, the noise in two hanning tapered 
channelmaps may be correlated depending on their sep- 
aration. Two channelmaps separated by one channelmap 
are correlated because both contain a quarter of the flux 
from the channel map between them. Only channel maps 
separated by more than one channel are independent. This 
will be shown in the following three cases in which two 
hanning tapered channel maps at different separations will 
be added. 

1. Adding two adjacent hanning tapered channelmaps i 
and (i+1) gives a signal H( i ) + ( i+1 ) of 

H(i)+(i+i) = Hi + H i+ i 

= fiUi-! + \Ui + \u i+1 ) 
+ (\u t + \u l+l + \u l+2 ) 

= jUi-i + jUi + I Ui+i + \Ui + 2 
and the noise cr ( l i ) + ( i+1 ) m that map will be 

+(R + i) 2 + (R +2 ) 2 ]^ 

— rr u — ^20 4 h 

- 4 4 y^g 

= \[i\(T h = 1.83 a h 

2. Adding the hanning tapered channels i and (i+2) gives 

H(i)+(i+2) = Hi + H i+2 

= {\Ui^ + \Ui + \u i+1 ) 
+ {\U l+l + \U i+2 + \U i+ z) 

= \Ui-\ + \\Ji + 5^i+l 
+ \Ui+2 + jU i+ 3 

and the noise becomes 

h r (1'u, \2 I /l u\2 I / 1 _u \2 

(T ( 4 ) + ( J +2) = [ + (.2 a i ) + (2^+1) 

+(i^ 2 ) 2 + (K +3 ) 2 ]^ 

vTS u Via 4 h 
= — a =—7E a 

= y^cr' 1 = 1.53 (J 71 

3. Adding the hanning tapered channels i and (i+3) gives 

H(i)+(i+3) = Hi + H i+3 

= (\U i - 1 + ^U i + ^U i+1 ) 
+ (\U i+2 + |^+3 + \U i+i ) 
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with a resulting noise of 



7 W + (*+3) 



(K-i) 2 + (K) 2 + (R+i) 2 
KK +2 ) 2 + (K + 3) 2 + (K + 4) 2 ]* 

/12 _ ^/l2 4 



4 " 

V2a h 



4 V6 



1.41 cr n 

So, channelmaps i and (i+3) are independent. 

Because the noise is correlated, adding N adjacent hanning 
tapered channelmaps does not give an increase of the noise 

with a factor \HSf but with a factor \ N — § • -4= as is 



4 

shown below. First the total signal Hat is calculated. 



Channel 


U(_i 


U, U i+ i ■■ 


Ui + jV-2 


Uj+AT-1 




i 

i+1 
i+2 


1/4 


1/2 1/4 
1/4 1/2 

1/4 ■■ 








i+N-3 
i+N-2 
i+N-1 






■ 1/4 
1/2 
1/4 


1/4 
1/2 


1/4 + 






|U 4 U 4+ i ■• 


Ui + N-2 


|U, +JV -i 





and thus 



Hjv — jUj-i + |Ui + Uj+i + • • • + Ui+N-2 + 



-A" 



From this it follows that the noise ct^ is given by 



The average continuum map to be subtracted is then 
formed by 

< C >= |(C/o«; + Cfiigh) 

Since crc !o „, and (T^. h are independent it follows that the 
noise <7<c> in the finally averaged continuum map is given 

by 



o<c> 



= 



AN 2 



+ 



N 2 



AN 2 



After subtraction of the continuum the channelmaps only 
contain signal from the HI emission line. The signal in the 
channelmaps containing the line emission is now given by 



U = Hi- < C> 

= \Ui-i + \Ui + \U i+1 - < C> 

Because a < c> is independent from af in the velocity 
range which is not used to form the averaged continuum 
map, it can be written 



N 



[ ( l)2 + ( 3 ) 2 + (Ar _ 2) . 1 2 + ( 3 ) 2 + ( l )2] - 



= [AT- §p v«=[N-?\ 



4 Ji 



(A - |)B' 1 o h 



However, before the hanning tapered channelmaps are 
added to form a total HI map, the continuum must be 
subtracted. This operation introduces extra noise in the 
channelmaps which doesn't behave like a hanning tapered 
correlation. Here, it will be assumed that the average con- 
tinuum map is formed by averaging Ni line free channels 
at the low velocity end of the datacube and N 2 channels 
at the high velocity end which gives 



1 Nl 

Q ow = — H . 



and 



1 ^2 



Since all channels are hanning tapered the noise in these 
maps can be calculated according to 



J, \(~L„u \2 , (lu\2 , i\ u \2 , J2 15 
a i = [\4 a i-l> +(2 a i) +\4 a i+l> +(7 <C>\ 

= V(i+^ 2 ) ° u 

When adding N adjacent hanning tapered and continuum 
subtracted channelmaps containing the line emission, the 
signal Lat will be 

Ln = jUi-i + jUi + Ui+i + ■■ ■ 

• • • + Ui+N-2 + jUi + N-l + \Ui + N 

-N-<C> 

The noise a l N at each pixel in the final map can be 
derived analogous to the calculation of cr^ and is given by 

a l N = [(N - §) + N 2 N 2 ] h <j u 




\) + N 2 M 2 



B h 



and 
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Fig. B.l. Integrated HI maps of all galaxies in the cluster more inclined than 45 degrees and brighter than M(B)=- 
16.5. The angular resolution is 30x30 arcsec. Contour levels are at column densities of 0.5, 2.0, 3.5 and 5.0 xlO 21 
atoms cm~ 2 . Individual galaxies are four times enlarged. To avoid overlap, some galaxies are slightly displaced from 
their actual position. N3906, N3938, U6962 and U6930 are more face-on than 45 degrees. U6940 is fainter than the 
complete sample limit. The dashed circle indicates the adopted boundary of the cluster. 
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Fig. B.2. Same as figure 1, displayed here in grayscale. Note the interacting pairs N3769/1135+48, N3893/N3896 and 
U6962/U6973. The galaxies N3718, N3726, N4010, N4013, N4088 and N4138 are strongly warped. It is clear that the 
integrated column density strongly depends on inclination. Many of the more face-on galaxies show a depletion of the 
HI gas in their inner regions. 



Observing parameters for UGC 6399 



Length of observation 


(hours) 


1 x 12 


Dates of observation 




/: a ,,„nn 
oAug^U 


Field center, a(1950) 




1 l.ZU.JU 


5(1950) 




51:10:00 


Central frequency 


(MHz) 


141D.JD 


Vhei of central channel 


(kms 1 ) 


810 


Primary beam FWHM 


(arcmin) 


37.4 


Nr. of interferometers 






Baselines (min-max-incr) (m) 


9700 79 


Synthesized beam (ax<5)(arcsec) 


12.9x 16.6 


Bandwidth 


(MHz) 


2.5 


Number of channels 




127 


Channel separation 


(kms" 1 ) 


4.15 


Velocity resolution 


(kms" 1 ) 


8.29 


rms noise in one channel (K) 


5.54 


K-mJy conversion, 






equiv. of lmjy/beam 


(K) 


2.84 



Contour levels for U6399 



Results from WSRT data 



From continuum map: 
21-cm flux density: 

central point source (mJy) <1.3 (3<r) 

extended source (mJy) <2.5 (3<r) 
From global profile: 

Integrated Hi-flux (Jy kms 4 ) 10.5 ± 0.3 

Hel. systemic velocity (kms -1 ) 791.5 ± 0.6 

HI profile width, 20% (km s" 1 ) 188.1 ± 1.4 

50% (kms 4 ) 172.5 ±2.9 
From velocity field: 

Hel. systemic velocity (kms -1 ) 794.5 ± 1.3 

Dynamical center, a(1950) 1 1:20:35.8 

5(1950) 51:09:58 

From total HI map: 

Geometric center, a(1950) 11:20:35.7 

5(1950) 51:10:01 

Position angle (deg) 139 

Inclination angle (deg) 77 

Diameter of HI disk (arcmin) 3.4 



Channel maps: 

(7=3.55 (K) 

Raw continuum map: 

(7=1.35 (K) 
Cleaned continuum map: 

(7=1.19 (K) 
Position- Velocity diagrams: 

(7=3.12 (K) 

Velocity fields: 

794.5 ±nx 15 (kms^ 1 ) 
Residual velocity field: 

± nx5 (kms~') 

Integrated HI map: 
0.52, 1.04, 

1.56 (x 10 21 atoms cm 4 ) 



Note: This is the most isolated system 
in the cluster. The rotation curve does 
not show a flat part on the receding side 
and maybe just barely on the approach- 
ing side. 
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Channel maps at a resolution of 12" x 16" x 19 km s~'. 
Contour levels at -3, -1.5 (dashed), 1.5, 3, 4.5, ... x(7. 
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